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Abstract. It is shown that a band hmited function on a non-compact sym- 
metric space can be reconstructed in a stable way from some countable sets 
of values of its convolution with certain distributions of compact support. A 
reconstruction method in terms of frames is given which is a generalization of 
the classical result of DufSn-Schaeffcr about exponential frames on intervals. 
The second reconstruction method is given in terms of polyharmonic average 
splines. 



1. Introduction 

One of the most interesting properties of the so called band limited functions, 
i. e. functions whose Fourier transform has compact support, is that they are 
uniquely determined by their values on some countable sets of points and can be 
reconstructed from such values in a stable way. The sampling problem for band 
limited functions had attracted attention of many mathematicians [3] , [4] , [5] , [14] . 

The mathematical theory of reconstruction of band limited functions from dis- 
crete sets of samples was introduced to the world of signal analysis and information 
theory by Shannon [21]. After Shannon the concept of band limitedness and the 
Sampling Theorem became the theoretical foundation of the signal analysis. In the 
classical signal analysis it is assumed that a signal (= a function) propagates in 
a Euclidean space. The most common method to receive, to store and to submit 
signals is through sampling using not the entire signal but rather a discrete set of 
its amplitudes measured on a sufficiently dense set of points. 

In the present paper we consider non-compact symmetric spaces which include 
hyperbolic spaces. A sampling Theory on Riemannian manifolds and in Hilbert 
spaces was started by the author in [16]- [20]. It was shown in [16] that for properly 
defined band limited functions on manifolds the following sampling property holds 
true: they can be recovered from countable sets of their amplitudes measured on a 
sufficiently dense set of points. 

Recently A. Kempf [12], [13] explored our results [16], [18] to developed a theory 
which approaches quantization of space-time and information through the sampling 
ideas on manifolds. 

By using Harmonic Analysis on symmetric spaces we give further development of 
the theory of band limited functions on manifolds. It is shown that under certain 
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assumptions a band limited function is uniquely determined and can be recon- 
structed in a stable way from some countable sets of values of its convolution with 
a distribution. 

The problem of deconvolution is very natural in the context of Signal Analysis 
since convolution is used for dcnoising a signal with subsequent reconstruction of 
the band limited component. On the other hand, since it is impossible to perform 
an exact measurement of a signal at a point, convolution with a distribution can 
represent a measuring device. 

As an illustration of the main result it is shown in Section 7 that every band 
limited function is uniquely determined and can be reconstructed in a stable way 
from a set of averages of its "derivatives" A"/, n G N1J{0}, over a set of small 
spheres of a fixed radius whose centers form a sufficiently dense set, here A is 
the Laplace-Beltrami operator of X. As a particular case we obtain that every 
band limited function on a non-compact symmetric space of rank one is uniquely 
determined and can be reconstructed in a stable way from a set of samples of its 
"derivatives" A"/(a;i),n e N1J{0}, for a countable and sufficiently dense set of 
points {xi}. Another particular case is, that a band limited function / on X is 
uniquely defined and can be reconstructed in a stable way from a set of its average 
values over a set of small spheres whose centers form a sufficiently dense set. 

The main result of Section 7 is the Theorem 7.3 which shows that for all these 
results the density of centers of the spheres depends just of the band width co and 
independent of the radius of the spheres r and of the smoothness index n. The 
radius of the spheres r depends on oj and smoothness n but it is independent of the 
distance r between centers and can be relatively large compare to r. In other words, 
a stable reconstruction can still take place even if spheres intersect each other. 

Reconstruction of band limited fimctions in L-2{R'^) from their average values 
over sets of full measure were initiated in the paper [6]. Our approach to this 
problem is very different from the approach of [6] and most of our results are new 
even in the one-dimensional case. It seems that reconstruction from averages over 
sets of measure zero (other then points) was never considered for band limited 
functions on W^. 

Let G be a semi-simple Lie group with a finite center and K its maximal compact 
subgroup. We consider a non-compact symmetric space X = G/K. We say that a 
function / from L2{X) is w-band limited if its Helgason- Fourier transform f{X,b) 
is zero for < A, A >-^/^= ||A|| > w, where < .,. > is the Killing form. We discuss 
some properties of such functions in the Section 3. It is shown in particular, that 
for any w > and any ball U C X restrictions of all w-band limited functions to 
the ball U are dense in L2{U, dx). 

In Section 4 we prove uniqueness and stability theorems. In Section 5 a recon- 
struction algorithm in terms of frames is presented. This result is a generalization 
of the results of Duffin-Schaeffer [5] about exponential frames on intervals. Let us 
recall, that the main result of [5] states, that for so called uniformly dense sequences 
of scalars {xj} the exponentials {e"^'>} form a frame in L2{I) on a certain interval I 
of the real line. In other words, if the support of the Fourier transform of a function 
/ e L2{R) belongs to I, then / is uniquely determined and can be reconstructed 
in a stable way from its values on a countable set of points {xj}. 

In Section 6 we consider another way of reconstruction by using spline-like func- 
tions on X. It is a far going generalization of some results of Schoenberg [22]. 
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2. Harmonic Analysis on symmetric spaces 

A Riemannian symmetric space X is defined as G/K, where G is a connected 
non-compact semi-simple group Lie with Lie algebra with finite center and K its 
maximal compact subgroup. Their Lie algebras will be denoted respectively as 
g and k. The group G acts on X by left translations and it has the "origin" 
o — eK, where e is the identity in G. Every such G admits Iwasawa decomposition 
G = NAK, where nilpotent Lie group N and abelian group A have Lie algebras 
n and a respectively. The dimension of a is known as the rank of X. Letter M is 
usually used to denote the centralizer o{ A in K and letter B is commonly used for 
the factor B = K/M. 

Let a* be the real dual of a and W be the Weyl's group. The S will be the set 
of all bounded roots, and S+ will be the set of all positive bounded roots. The 
notation a"*" has the following meaning 

a+ = {/i e a\a{h) > 0, a G E+} 

and is known as positive Weyl's chamber. Let p £ a* is defined in a way that 2p is 
the sum of all positive bounded roots. The Killing form <, > on g defines a metric 
on a. By duality it defines a scalar product on a*. The a^ is the set of A G a*, 
whose dual belongs to a+. According to Iwasawa decomposition for every g G G 
there exists a unique A{g) e a such that 

g = n exp A{g)k, k G K,n G N, 

where exp : a — > A is the exponential map of the Lie algebra a to Lie group A. On 
the direct product X x B we introduce function with values in a using the formula 

(2.1) A{x,b) = A{u-^g) 

where x = gK, g G G,b = uM, u G K. 

For every / G Cq^{X) the Helgason- Fourier transform is defined by the formula 

f{X,b)= f /(x)e(-'^+'')^(^'''»dx, 
Jx 

where X G a* ,b G B = K/M, and dx is a G-invariant measure on X. This integral 
can also be expressed as an integral over group G. Namely, if 6 = uM, u G K, then 

(2.2) /(A,&)= / /(x)e(-^^+^)^(""'^))d5. 

JG 

The following inversion formula holds true 

fix) = w-^ I /(A, 6)e(*^+'')(^(^'''» \c{X)\-^dXdb, 

where w is the order of the Weyl's group and c(A) is the Harish-Chandra's function, 
dX is the Euclidean measure on a* and db is the normalized if-invariant measure on 
B. This transform can be extended to an isomorphism between spaces L2{X,dx) 
and L2(a+ x B, \c{X)\~'^ dXdb) and the Plancherel formula holds true 

11/11= (/ \f{XMM>^)\-^dXdb 

\JaL''_^x.B 

An analog of the Paley- Wiener Theorem is known which says in particular that 
a Helgason-Fourier transform of a compactly supported distribution is a function 
which is analytic in A. 
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To introduce convolutions on X we will need the notion of the spherical Fourier 
transform on the group G. 

For a A e a* a zonal spherical function (px on the group G is introduced by the 
Harish- Chandra's formula 

(2.3) </?a(.9) = / e'''^+P^^^^^3))^k 

Jk 

If / is a smooth bi-invariant function on G with compact support its spherical 
Fourier transform is a function on a* which is defined by the formula 



/(A) = / f{g)^-x{g)dg. 
Jg 

The inversion formula is 

J a" 

The corresponding Plancherel formula has the form 



G 



1/2 / ^ 

2, 



1/2 



1/(5)^5 = / IfWl'dX 



+ 



If / is a function on the symmetric space X and V' is a K-bi-invariant function 
on G their convolution is a function on X which is defined by the formula 

/ * ^(.9 • o) = / f{gh-^ ■ o)^{h)dh, g e G. 



JQ 

By using duality arguments the last definition can be extended from functions to 
distributions. It is known, that 

(2.4) /7^(A,6) = /(A,6)V^(A). 

Denote by Tx{M) the tangent space of M at a point x G M and let exp^ ■ 
Tx{M) ^ M be the exponential geodesic map i. c. expj.{u) = 7(1), m € Tx{M) 
where 'y{t) is the geodesic starting at x with the initial vector u : 7(0) = x, ^^^^ = u- 
In what follows we assume that local coordinates are defined by exp. 

By using a uniformly bounded partition of unity {ifv} subordinate to a cover of 
X of finite multiplicity 

X = \jB(x,,r), 

where B[xv,r) is a metric ball x^ G X of radius r we introduce Sobolev space 
H''{X),a > 0, as the completion of C^{X) with respect to the norm 

(2-5) II/IIh^W= (Ell'^-/llH'(B(..,r))) • 

The usual embedding Theorems for the spaces H'^{X) hold true. 

The Killing form on G induces an inner product on tangent spaces of X. Using 
this inner product it is possible to construct G-invariant Riemannian structure on 
X. The Laplace-Beltrami operator of this Riemannian structure is denoted as A. 

If the space X has rank one, then in the polar geodesic coordinate system 
{r,0i, Od-i) on X at every point x G X it has the form [9] 

A = + -^r\^ Or + As, 

b (r) dr 
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where As is the Laplace-Behrami operator on the sphere S{x, r) of the induced 
Riemannian structure on S{x, r) and S{r) is the surface area of a sphere of radius 
r which depends just on r and is given by the formula 

(2.6) S{r) = Q.i2-'i c-P-'ishP{cr)sh'i{2cr), 

where d = dimX = p -\- q + l,c = {2p + Sg)"^/^, p and q depend on X and 

2^-^/2 



V{d/2) 



is the surface area of the unit sphere in d-dimensional Euchdean space. 

In particular, if a function / is zonal, i.e. depends just on distance r from the 
origin o, we have 

(2.7) M{r) = 4:rAis{rf^'^ 



S{r) dr \ dr 

Since for a smooth function / with compact support the function A/ can be 
expressed as a convolution of / with the distribution A6e, where the Dirac measure 
Se is supported at the identity e of G, the formula (2.4) implies 



(2.8) A/(A,6) = -(II Af + \\pf)f{X,b),f€ C^iX), 

where ||Ap =< A, A >, ||p||2 =< p,p >,<,> is the Killing form. 

It is known that for a general X the operator (—A) is a self-adjoint positive 
definite operator in the corresponding space L2{X, dx), where dx is the G-invariant 
measure. The regularity Theorem for the Laplace-Beltrami operator A states that 
domains of the powers (— A)'^/^ coincide with the Sobolcv spaces H'^{X) and the 
norm (2.5) is equivalent to the graph norm ||/|| + ||(— A)'^/2/||. 



We consider a ball B{o,r/A) in the invariant metric on X. Now we choose 
such elements g„ & G that the family of balls B{xv,r/4),Xv = ■ o, has the 
following maximal property: there is no ball in X of radius r/4 which would have 
empty intersection with every ball from this family. Thcin the; balls of double radiiis 
B{x^,r/2) would form a cover of X. Of com-se, the balls B{xv,r) will also form a 
cover of X. Let us estimate the multiplicity of this cover. 

Note, that the Riemannian volume B{p) of a ball of radius p in X is independent 
of its center and is given by the formula 

B(p)= r s{t)dt, 

Jo 

where the surface area S{t) of any sphere of radius t is given by the formula (2.4). 

Every ball from the family {i?(.T,y, r)}, that has non-empty intersection with a 
particular ball B{xj, r) is contained in the ball B{xj, 3r). Since any two balls from 
the family {B{x„, ''/4)} are disjoint, it gives the following estimate for the index of 
multiplicity of the cover {B{xv,r)}: 

^(3r) _ CS{t)dt 
^ ■ ' Bir/4) J^/'s{t)dt 



6 



ISAAC PESENSON 



It is clear that for all sufficiently small r > this fraction is bounded. In what 
follows we we will use the notation 

N = sup , ,\ . 
o<r<i B{r/4) 

So, we proved the following Lemma. 

Lemma 2.1. For any sufficiently small r > there exists a set of points {x/j,} from 
X such that 

1) halls B(Xfj,,r/4) are disjoint, 

2) balls B{x^,r/2) form a cover of X, 

3) multiplicity of the cover by halls B{Xfj,,r) is not greater N. 

We will use notation Z({a;^}, r, TV) for any set of points {a;^} G X which satisfies 
the properties 1)- 3) from the last Lemma and we will call such set a metric (r, A^)- 
lattice of X. 

The following results can be found in [16] for any homogeneous manifold X. 

Theorem 2.2. For any k > d/2 there exist constants C = C{X,k,N) > 
0, ro{X, fc, A^), such that for any < r < ro and any (r, N)-lattice Z = Z({a;^}, r, A'') 
the following inequality holds true 



(2.10) 11/11 <c|r'^/2 |/(x,)pj +r'=||A'=/V|l|,fc>rf/2. 

and there exists a constant Ci = C\{X,k,N) such that 

1/2 

E <C4f\\Hu^x),f&H\X). 

3. Band limited functions 

Definition 1. We will say that / G L2{X,dx) belongs to the class B^{X) if its 
Helgason- Fourier transform has compact support in the sense that f{X,h) = for 
||A|| > CO. Such functions will be called w-band limited. 

We have the following important result, which is a specification of some results 

in [15]- [18]. 

Theorem 3.1. A function f belongs to B^iX) if and only if the following Bernstein 
inequality holds true for all a > 0, 

(3.1) ||A'^/||<(a;2 + ||p|p)-||/||. 

Proof. By using the Plancherel formula and (2.8) we obtain that for every cj- band 
limited function 

v<T f\\2 _ I I ^||\||2 I ||„M2^(7| 

||A||<w Jb 



\\xr+\\prr\mwmrdxdb< 

+ \\pfr [ [ l/(A, b)f\c{X)\-'dXdb = {iv' + \\pfy 
Ja* Jb 



\\2\an fn2 
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Conversely, if / satisfies (3.f), then for any e > and any a > we have 
/ / \f{X,b)\M>^)\-^dXdb< 

J\\X\\>U+£ Jb 

I I (iiAf + hrr'^m? + \\pff''\f{xMM^)\-^dxdh < 

J\\\\\>u+e Jb 



(3.2) 



\pf 



la 



|2 



It means, that for any £ > the function /(A, 6) is zero on {A : ||A|| > w + e} x B. 

The statement is proved. □ 

Now we are going to prove the following "density" result. 

Theorem 3.2. For every w > and every ball U C X restrictions to U of all 
functions from Bjj{X) are dense in the space L2{U,dx). 

Proof. Indeed, assume that ip € L^iUjdx) is a function which is orthogonal to 
all restrictions to U of all functions from B^{X). We extend by zero outside 
of U . By the Paley- Wiener Theorem the Helgason-Fourier transform il}{X, h) is 
holomorphic in A and at the same time should be orthogonal to all functions from 
L2 (B(0,w) X B; \c{X)\-'^dXdb), where 

S(0,a;) = {AGa* : ||A|| < w}. 

It implies that 'ij) is zero. Consequently, the function ^ is zero. It proves the 
Theorem. □ 

4. Uniqueness and stability 

We say that set of points M = {xj] is a uniqueness set for B^^iX), if every 
/ G B^{X) is uniquely determined by its values on M. 

Theorem 4.1. If a set M = {xj} is a uniqueness set for the space B^{X), then 
for any bi-invariant distribution of compact support ^ every function f G B^{X) is 
uniquely determined by the set of values f * (j){xj ) . 

Proof. If f & Bi^{X) then because 

f*<j) = f(j) 

the function f *(j} also belongs to Bi^{X) and by assumption is uniquely determined 
by its values / * (pi^j)- But by the Paley- Wiener Theorem the zero set of the 
function (j) has measure zero. Thus, the equality 

/i * ^ - /a * = 0(7i - fi) = 0, 
implies that /i = /2 as L2-functions. The statement is proved. □ 

For any uniqueness set M and any w > the notation Z2 (-^) will be used for 
a linear subspace of all sequences {vj} in I2 for which there exists a function / in 
Bi^{X) such that 

/ * (p{xj ) = Vj, Xj S M. 

In general 12{M) ^ I2. 



8 



ISAAC PESENSON 



Definition 2. A linear reconstruction method R from a uniqueness set M is a 
linear operator 

R : Z^(M) ^ B^{X) 

such that 

The reconstruction method is said to be stable, if it is continuous in topologies 
induced respectively by I2 and L2{X). 

The Theorem 4.1 does not guarantee stability but the next one does. 

Theorem 4.2. There exists a constant c = c{X,N) > such that for any u! > 0, 

any {r, N)-lattice Z{{x^},r,N) with 

Q<r<c{oj^ + \\pf)-^'\ 

and every bi-invariant distribution (j) with compact support, whose Helg as on- Fourier 
transform </> does not have zeros on B(0,w), every function f G is uniquely 

determined and reconstruction method from the set of samples f * (t>{xj) is stable. 

Proof. The formulas (2.10) and (3.1) imply, that there exist constants C, c such 
that for any w > 0, any Z = Z{{x^}, r, N) with 

< r < c(u;2 ^ ||p|n-V2 

and every / G B„(X) 

1/2 



(4.1) 



Indeed, by the Theorem 2.2 and the Theorem 3.1 we have 



(4.2) 



<C { 



1/2 

E \f{x,)\^\ +r\u:- + \\prf/-\ 



where k > d/2. If we choose 

0<r <c{uj^ + \\p\\Y^/\c = C-\ 

then the inequality (4.2) gives the inequality (4.1) which implies uniqueness. Next, 
by applying the Planchcrcl Theorem and the inequality (4.1) to the function f * (p 
we obtain the inequality 

11/11 = ll7ll = lir'/?ll<c^i||/^ll< 



1/2 



(4.3) 
where 



\XjeZ{xj,r,N) 



inf (j)(X) 



This inequality (4.3) implies stability of the reconstruction method from the samples 
/ * (l>{xj). The Theorem is proved. □ 
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5. Reconstruction in terms of frames 

In this section wc are going to present a way of reconstruction of a function 
/ G Bu:{X) from a set of samples of its convolution f * (p in terms of frames in 
Hilbert spaces [1], [5]. 

If is a Dirac distribution at a point Xj & X then according to the inversion 
formula for the Helgason-Fourier transform we have 

<<5x,,/> = w-^ [ /(A,6)e(^^+'')(^(^-*»|c(A)|-2rfAd6. 

Ja'xB 

It implies that if / G L2{X) then the action on f{X,b) of the Helgason-Fourier 
transform of is given by the formula 

(5.1) f{X,b) ^ (C;,/) = [ e(^^+'')(^(--''))/(A,6)|c(A)|-2dAd&. 

Theorem 5.1. There exists a constant c = c{X, N) such that for any given u) > 0, 
for every (r, N)-lattice Z = r, N) with 

< r < c(u;2 ^ llpf )-i/2, 

the following statements hold true. 

1) The set of functions {5xj} forms a frame in the space 

(5.2) La (S(0,a;) x B- \c{\)\-'^ dXdh) 

and there exists a frame {Oj} in the space Bi^{X) such that every co-band limited 
function f G B^{X) can be reconstructed from a set of samples {5xj{f)} by using 
the formula 

(5.3) f=Y,6x,{f)Q,. 

2) If the Helgason-Fourier transform of a com,pactly supported distribution (f> 
does not have zeros on the ball B{Q,u)), then the Helgason-Fourier transforms of 
the distributions 

f f*(j){Xj) 

form a frame in the space L2 (/B(0,cj) x B; |c(A)|~^(iA(i&) and there exists a fram,e 
{^j} in the space B^^X) such, that every f G Bi^{X) can be reconstructed from 
the samples of the convolution f * (j) by using the formula 



(5.4) /= ^ ^,^.(/*^)$^.. 

Proof. An application of the Plancherel formula along with the inequality (4.3) and 
the second inequality in the Theorem 2.2 give that for any / e Bi^{X) 

(5.5) Ci||7||a.< f E I<^'^'/>I') <^^2||7||a., 

where < .,. > is the scalar product in the space A2 = L2 (B(0,ti;) x B\\c{\)\~'^d\db^. 
The first statement of the Theorem is just another interpretation of the inequalities 
(5.5). 
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We consider the so called frame operator 

j 

It is known that the operator F is invertible and the formula 

(5.6) 6^ = ^-1^;:. 

gives a dual frame Qj in L2 [B{0,oj) x B; \c{X)\~'^dXdb) . A reconstruction formula 
of a function / can be written in terms of the dual frame as 

(5.7) / = 5Z<e;,/>^ =5]<C;.,/>g;, 

j j 

where inner product is taken in the space L2 (S(0,w) x B; \c{X)\~'^dXdb) . 

Taking the Helgason- Fourier transform of both sides of the last formula we obtain 
the formula (5.3) of our Theorem. 

The second statement of the Theorem is a consequence of the first one, of the 
Plancherel formula and of our assumption that (j) does not have zeros in B{0,uj). It 
is clear, that in the formula (5.4) every $j is the inverse Helgason-Fourier transform 
of the function □ 

In the classical case when X is the one-dimensional Euclidean space we have 

(5.8) C;.(A)=e-^\ 

In this situation the first statement of the Theorem means that the complex expo- 
nentials e'^^^^,Xj e Z{{x^},r, N) form a frame in the space L2{[—uj,oj]). 

Note that in the case of a uniform point- wise sampling in the space L2{R) this 
result gives the classical sampling formula 

m = E f(^-^) ^ft-^Io) ' " = ^^z^' ^ < ^' 

with a certain oversampling. 



6. Reconstruction by using polyharmonic splines on X 

To present our second method of reconstruction we consider the following op- 
timization problem. Although the results of this section are similar to the corre- 
sponding results from our paper [17], the presence of the Helgason-Fourier transform 
allows to make all our constructions much more explicit. 

For a given {r, N) -lattice Z = Z{{Xfi},r,N) find a function in the Sobolev 
space H'^'^{X) for which L'^iXfj) = Sv^, € Z, and which minimizes the functional 
u — )■ ||A'^u||. Here the d^f^ is the Kronecker delta. 

Theorem 6.1. For a given (r, N)-lattice Z{{Xfj_},r, N) the following statements 
are true. 

1 ) The above optimization problem does have a unique solution. 

2) For every there exists a sequence a = {aj>} S I2 such that 

(6.1) Lt{\b)=Y,o^3JA\b), 

Xjez 
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where 6j is the Helgason-Fourier transform of the distribution Sj and is given by 
the formula 

(6.2) 4, = e(-^^+'')(^(^^'''»,Xj- G Z. 



Proof. The Theorem 2.2 imphes that for a fixed set of values {vj} G h the minimum 
of the functional 

u ||A''u|| 

is the same as the minimum of the functional 

1/2 



2 

'J I 



L 



with a set of constrains u(xj) = Vj, Xj G Z{{Xf^}, r, N). 

Since the last functional is equivalent to the Sobolev norm it allows to perform 
the following procedure. 

For the given sequence Vj = 5j,i, where v is fixed natural number consider a func- 
tion / from if^'^(X) such that f{xj) = vj. Let Pf denote the orthogonal projection 
of this function / (in the Hilbert space H'^'^{X) with natural inner product) on the 
subspace of functions vanishing on Z{{Xf^}, r, N). Then the function = f — Pf 
will be the unique solution of the above minimization problem. 

The condition that G H'^'^[X) is a solution to the minimization problem 
implies, that A'^Lj should be orthogonal to all functions of the form A'^/i, where 
h G H'^^{X) and has the property h{xij) = for all ji. This leads to a differential 
equation 

{is!" L'D'EFhdx = '^aj^Jiixj), 

j 

for an I2- sequence {aj^^}. Thus, in the sense of distributions 

(6.3) A2'=L^ = ^a,.A^. 

3 

Taking the Hclgason-Fouricr transform of both sides of (6.3) in the sense of 
distributions we obtain the equation (6.1). 

The Theorem is proved. □ 

We will need the following Lemma. 

Lemma 6.2. If for some f G H'^'^{X),a,a > 0, 

(6.4) 11/11 <a||A-/||, 
then for the same /, a, a and all s >0,m = 2', I = 0,1, 

(6.5) IIAVII < a-IIA^^+VII, 

iffe fl'2(m<T+s)(^)_ 

Proof. In what follows wc use the notation d^i{X) for the measure |c(A)|~^dA. The 
Plancherel Theorem allows to write our assumption (6.4) in the form 

ll/f=/ |/(A,6)|2d^(A)d6< 

Ja+xB 
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a' / (||A||2 + ||pf)2-|/(A,6)pdMA)d6< 

Ja+xB 

[ [ i\\Xf + \\pfr'^\f{\,b)\-'d,,{X)db+ 

a' f I {\\\r + \\Prf''\h\b)\'dn{X)db. 

•''(l|A|P + ||p|P)>a-V- Jb 



^P+IIpP) 
From this we obtain 



\f{\h)\^ - a^dlAf + \\p\\y'^\f{\b)f) dn{X)db 



0< / / {\f(X,b)\' -a'(\\Xr + \\prY^\f(X,bW]dti(X)db< 

•^(l|A|P+||p|P)<o-i/- Jb 

[ [ (a'iWXf + \\pfr\f{X,b)\' - |/(A,6)|2) rfM(A)rf6. 

Multiplication of this inequality by a^(||A||^ + will only improve the 

existing inequality and then using the Plancherel Theorem once again we will obtain 

11/11 <a||A<^/||<a^||A2VII- 
It is now clear that using induction we can prove 

11/11 <a-||A™VI|,m = 2',ZeN. 
But then, using the same arguments we have for any s > 

0< / / (a'^{\\Xf + \\pff^'^\f{X,b)\'- 

a2(m+.)(p||2^ ||p||2)2(™+.)a|;(^^^)|2)^^(^)^^ < 

f f (a2(m+.)(||A||2^||^||2)2(m+sV|^'(^^^)|2_ 

•^(l|A|P + ||pP)>a-i/- Jb 

a'-m\xr + \\prr^\.f{x,b)\'')dp{x)db, 

that gives the desired inequality (6.5) if t = as. □ 



Lemma 6.3. There exists a constant C = C{X,N), and for any co > there 
exists a ro{X,N,uj) > such that for any < r < ro{X,N,co), any {r, N) -lattice 
Z = Z{{xn},r,N), the following inequality holds true 



\J2fi^j)Lt'^'-f\\HsiX)< 



(6.6) {Cr'{u;' + \\pf))'''{co' + \\p\\r\\f\\, 

for any s > 0,1 = 0,1, and any f G Bu{X). 
Moreover, if s > d/2 + k then 

||^/(x,)if+«-/||c,.(^)< 
j 

(6.7) {Cr'iu;' + Wpf))''" {u' + \\pfr ll/IM = 0, 1, .•• 

where C^{X) the space of k continuously differentiable bounded functions on X. 
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Proof. First we will show that there exist constants C = C{X, N), and ro{X, N) > 
such that for any < r < rQ{X,N), any (r, A'')-lattice Z = Z{{xn},r,N), the 
following inequality holds true 

(6.8) \\Y;^f{xj)Lf''+^ - < 2{Cr'Y\\^''''^'fl 

3 

for any s > 0,1 = 0,1, and any / G H'^{X). 

Indeed, by the Theorem 2.2, since for = d the function Ylj f{^j)L] '^'^^ inter- 
polates / we have 

\\Y,f{x,)L^;'+^ - /II < cr^'w^^Y.f^'^i)^"^'^' - 

3 3 

By the Lemma 6.2 we obtain 

||A^(^/(x,)L,^''^+^-/)il < {Cr^'fW+\Y.fi^^)Ll'^' - 

3 3 

Using the minimization property we obtain (6.7). 

li s > d/2-\-k then an application of (6.7) and the Sobolev embedding Theorem 
gives 

(6.9) II ^ /(a;,)Lf - /||c,.(^) < {Cr^'f W^+^fll = 0, 1, 

3 

where C^{X) the space of k continuously differentiable bounded functions on X. 
The inequalities (6.7) and (6.8) along with the inequality 

||A'^/||<(a;2 + ||p||y 11/11, /eB.(X), 

imply the Lemma. □ 

Inequalities (6.6) and (6.7) show that a function / € i3^(X) can be reconstructed 
from its values f{xj) as a limit of interpolating splines 



5:/*<^(a.,)Lf'^+^ 



in Sobolev or uniform norms when I goes to infinity. By using the Planchercl 
formula for the Hclgason-Fourier transform wc obtain the following reconstruction 
algorithm in which we assume that a reciprocal of the Hclgason-Fourier transform 
of a distribution cj) is defined almost everywhere and bounded. 

Theorem 6.4. There exists a constant c = c{X,N), so that for any co > 0, any 

{r,N)-lattice Z = Z({a;^}, r, iV) with 

r<(c(..2 + |bf))-^/^ 

every function f G Bi^{X) is uniquely determined by the values of f*(j) at the points 
{xj} and can be reconstructed as a limit of 



^/*0(^,■)Af'*+^^^oo, 
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in Sobolev and uniform norms, where 

£?d+s 



A. 



2'd+s _ j 



7. Example: Convolution with spherical average distribution 

In this section we assume that the symmetric space X = G/H has rank one, 
which means that the algebra Lie a of the abeHan component A in the Iwasawa 
decomposition G = NAK has dimension one. For a point y G X, the S{y, r), r > 0, 
will denote the sphere with center y and of radius r. For a smooth function / with 
compact support we define the spherical average 

{M^f){y) = S{T)-' f f{z)ds{z),r>0, 

Jsiy.r) 

where ds is the measure on S{y, r) and S{t) is its volume, i.e. the surface integral 
of ds over S{t) which is independent on the center. 

It is clear, that M'^ f is a convolution of / with the distribution 

m^{f) = S-Hr) f f{x)ds{x),T > 0,/ e G^{X). 

Js{o,t) 

Since X is a symmetric space of rank one, the group K is transitive on every 
S{o,t) and this fact allows to reduce integration over the sphere S{o,t) to an 
integral over group K. Because the invariant measure dk on K is normalized, it 
results in the following change of variables formula 

(7.1) (MV)(5-o)= / f{9kzr)dk, 

Jk 

for any Zr G S{o,T),g G G. 

Next, by using the following formula 

/ f{x)dx= [ f{g-o)dgJ &G^{X) 
Jx Jg 

along with Minkowski inequality and the formula (7.1) we obtain 

1 /2 

||M-(/)|| <y^J^ \f{gkzr)\^dkdg^ < 

1 /2 

(7.2) Lilc}^^^"^^^'^^] dk=\\f\\,fec^{x), 

where we have used the fact that the invariant measure dk is normalized on K. 
Since the set C^{X) is dense in L2{X) it implies the following. 

Lemma 7.1. Operator f — > M'^{f) has continuous extension from C^{X) to 
L2{X) which will be also denoted as M'^{f) and for this extension 

\\M^f)\\ < 11/11 



for all f G L2{X),T> 0. 
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Since the operator M'^ is a convolution with the distribution rri^ , we have 

M^) = m-f. 

Our goal is to consider even more general operator. We choose an n G NU{0}, 
and consider the operator 

/^M-((-A)"/). 

It is clear, that this operator is a convolution with the distribution (— A)"m'^. 

The nearest goal is to find the Fourier transform {—A)""m'^ of the distribution 
(— A)"m^. Namely, we prove the following. 

Lemma 7.2. The Helgason-Fourier transform of the distribution (— A)"m'^ is 
given by the formula 

(7.3) (-A)^- = {\\Xf + llpf )>A(expry), 

where V is a vector from algebra Lie a and (p\ is the zonal spherical function i.e. 



{g)= [ e^*^+P^^^^''^^Uk,geG. 



Proof. Since the operator (— A)",n S NU{0}, can be represented as a convolution 
with the distribution A"^e) where 5^ has support at the identity e G G, it is clear, 
that 

(-A)^- = (||Af + 

But 

7^(A,6) =m^(e('^+^)^(^''')) = 

(7.4) S-^{t) [ e(*^+^)^(^'*)ds(x),T>0. 

JS{o.t) 

By the formula (7.1) and by the Harish-Chandra's formula (2.3) the last integral 

is 

S-\t) [ e^'^+P^^^^'''Us{x) = [ e(-'^+'')^('="'s)rffc = (^A(5), 
Js{o,t) Jk 

where x = gK, b = kM e B = K/M. 

Because the symmetric space X = G/K, has rank one, there exists a unique 

vector V that belongs to the Lie algebra a, such that the curve exp(— t^) • o is a 

geodesic in X and the one-parameter group exp(— rV) moves the origin a £ X to 

a point Zr in (7.1). It gives the formula 

m^{X,b) = fh^{X) = </?A(exp— rV). 

The Lemma is proved. □ 

The following result is a specification of the Theorem 5.1. 

Theorem 7.3. There exists a constant c = c{X, N) such that for any given u> > 
0,n e N|J{0}, for every {r, N)-lattice Z({x^},r,7V) with 

< r < c(u;2 + ||pf )-V2 

and any 

(7.5) r<(a;2 + ||pf)-("+^^/\ 
the following hold true. 
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1) There exists a frame € B^{X) such that for any f e Baj{X) the following 
reconstruction formula holds true 

(7.6) / = ^5.,((-A)"m. 

j 

which means that f can be reconstructed from averages of (—A)"/ over spheres of 

radius t with centers at Xj S Z({x^}, r, N). 

2) Every f € Bu{X) is a limit (in Sobolev and uniform norms) of the functions 

j 

when I goes to infinity and 

7-2'd+s 



(||AP + ||p||2)n^;,(exprF)' 
where i^'''"'"* is given by the formula (6.1) 

Proof. First we are going to show that the following estimate holds true 
KIIAf + ||pf)>A(exp(-Ty))-l|< 

(7.7) min{2(||Af + ||pf)";r2(||Af + ||pf )"+!}. 

Indeed, since tp^ is a spherical function, the function 

m?+\\pfTvx{eM-rV)) = ^x{T) 

is zonal and ^^\{t) can be calculated according to the formula (2.5). By using 

this fact along with the fact that the spherical function ^\ is an eigenfunction the 
Laplace-Beltrami operator (—A) with the eigenvalue (||A|p + ||/9|p)"+^ we obtain 

1-$a(t) = 



By using the following change of variables: u = sr, 7 = str, and the inequality 
Ifxig)] < Ij we obtain the estimate (7.7). 

According to inequalities (7.5) and (7.7) the Helgason-Fourier transform of the 
distribution (— A)"mr does not have zeros on the interval [— a;,cj]. Consequently, 
the first statement of the Theorem is a consequence of the Theorem 5.1. Note, 
that since the Helgason-Fourier transform of the distribution (— A)"TOt- is given 
by the formula (7.3) the Theorem 5.1 gives that every function $j is the inverse 
Helgason-Fourier transform of the function 

^ 



|A|P + ||pP)"'^A(expry)' 



where functions form a frame which is dual to the frame 6^^ in the space 
L2 ([-w,w] X B;\c{X)\-^dXdb) . 

The second part of the Theorem is a consequence of the Theorem 6.4. □ 
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As a particular case with t — we have the so called "derivative" sampling, 
which means that a function can be reconstructed from the values of A"/ as 

3 

for corresponding frame 

Another particular case we obtain when n = 0. In this situation the formula 
(6.6) tells that a function can be reconstructed from its averages over spheres of 
radius r with centers at Xj. 

The case t = 0, n = 0, corresponds to a point-wise sampling. 
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